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ABSTRACT 



OO 

^ ' Supermassive black hole binary (SMBHB) is expected with the hierarchical 

?-j ■ galaxy formation model. Currently, physics processes dominating the evolution 

^ ■ of a SMBHB are unclear. An interesting question is whether we could observa- 

tionally determine the evolution of SMBHB and give constraints on the physical 
processes. Jet precession have been observed in many AGNs and generally at- 
tributed to disk precession. In this paper we calculate the time variation of jet 



^ ■ precession and conclude that jet precession is accelerated in SMBHB systems 

but decelerated in others. The acceleration of jet precession dPpj./dt is related to 

dP P 

lO ■ jet precession timescale Ppr and SMBHB evolution timescale Ta, — ^ ~ —A—. 

■-■ , ^ at 

^ ■ Our calculations based on the models for jet precession and SMBHB evolution 

O ■ show that dPpr/dt can be as high as about —1.0 with a typical value —0.2 and 

can be easily detected. We discussed the differential jet precession for NGC1275 
^ ■ observed in the literature. If the observed rapid acceleration of jet precession 

■ is true, the jet precession is due to the orbital motion of an unbound SMBHB 

with mass ratio q ~ 0.76. When jets precessed from the ancient bubbles to the 
currently active jets, the separation of SMBHB decrease from about 1.46 Kpc 
to 0.80 Kpc with an averaged decreasing velocity da/dt ~ —1.54 x lO^cm/s and 
evolution timescale ~ 7.5 x 10^ yr. However, if we assume a steady jet preces- 
sion for many cycles, the observations implies a hard SMBHB with mass ratio 
q ^ 0.21 and separation a ^ 0.29 pc. 



Subject headings: accretion, accretion disks — galaxies: formation — galaxies: 
interactions — galaxies: individual (NGC1275(3C84)) — galaxies: jets — grav- 
itational waves 
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Introduction 



In the hierarchical galaxy formation models of cold dark matter (CDM) cosmology, 
present-day galaxies are the products of successive mergers. Recent observations show that 
almost all galaxies harbor at center a supermassive black hole (SM BH) of mass tightly 
corre l ating with both the mass and the velocity dispe r sions of the bulge (IFerrarese &: Merritt 



2000l : ICebhardt et al II2OO0I : iMagorrian et al. lllQQsl : IXremaine et al Il2002h . During galaxy 



interacting and merging, the g as at galactic plane is driv en toward central SMBH, triggering 
the activity of active galaxy (jWilson fc Colbert I Il995l ) and black hole accretion. SMBHs 
in galactic nuclei likely increase mainly through matter accretion. In this scenario, galaxy 
interacting and merging is expected to trigger the formations of an unbound binary active 
galactic nuclei (AGNs). In galaxy mergers, two galaxies and the SMBHs at center initially 
lose their orbital angular momentum owing to galactic dynamic friction and form a bound 
supermassive black hole binary (SMBHB) at a separation of an ~ 10 pc, when the SMBHB 
bind energy becomes dominated. The dynamics friction is very efficient because of the traps 
of stellar objects around each black hole and the evolution timescale of SMBHB is order of the 
local dynamic timescale, depending on the stellar velocity of two galaxies. The evolution of 
a bound SMBHB is dominated by the dynamic friction but the evolution timescale depends 
on the inner surface brightness profiles of galaxies. A SMBHB becomes hard at a separation 
Oh ~ 0.1 — 1 pc, when the loss of the orbital angular mome ntum is dominated by three-body 



interactions be tween SMBHB and the stars passing by (IBegelman et al. I Il980l : iQuinlan 



19961 : IYu II2OO2I ). When SMBHB becomes hard but orbital angular momentum loss because 
of gravitational wave radiation is unimportant, SMBHB may stall at a ~ Oh on a timescale 
longer than the Hubble time. However, observations of nearby galaxies suggest that most 
SMBHBs should have passed throug h the hard phase and become coalesced quickly. This 
is the so-called final parsec problem (IMerritt fc Milosavljevic II2005I ). To solve the problem, 
theoretically several processes w ith large uncertainties have been suggested in the literature 
( iMerritt fc Milosavlievic 1120051) and the hydrodynamics interaction with gas disk may play 



the maior role flGould fc Rix II2OOOI : iLiu et al. II2003I : iLiu II2004I : lArmitaee fc Nataraian 112005 



Escala et al. \\2004 ). An important question is whether we could detect SMBHB at galactic 
center, determine its evolution, and give observational constraints on the formation and 
evolution of SMBHB. 

Although unbound AGN binary systems with sep a.ration of order Kpc or larger have 



been imaged in interacting and merging galaxies (e.g. iKomossa et al. I l2003l : iBallo et al 
2004 : [Rodriguez et al. II2006I ). no hard or bound SMBHB has been directly detected. Close 
SMBHBs or binary coalescence have been introduced in e xplaining the observat i ons of many 



AGNs, for example, periodic optical and radi o outbursts (ISillanpaa et al. Ill988l : iKatz Ill997 



Liu et al. 1995 



19971 . l2006l : lLiu fc Wu 1120021 ). periodic variation of Very Long Baseline Inter- 
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ferometer (VLBI) jet position angle (IStirling et al. II2003I : ISudou et al. 1120031). th e interrup- 
tion of jet formation in double-doubl e radio galaxies (DDRGsl JLiu et al. 2003 ). X-shaped 



radio feature in winged radio sources (ILiu 



2004 



Merritt fc Ekers 1120021) . and S- or Z-shaped 



morphological symmetry of radio jets (jSegelman et al. Il980 . BBR ). A recent review on the 



Komossa 



J2006h . The periodic outbursts 



observational evidences for SMBHB was given by 
may be due to the periodic interaction of a SMBHB and a standard accretion disk or an 
ADAF, while the periodic variation of VLBI jet position angle is because of binary orbital 
motion. However, when the disk mass inside the binary orbit is less than the mass of the 
secondary which is at a radius order of 10^ times the Schwarzschild radius of the primary 
black hole, the interaction between the accretion disk and the SMBHB will real ign the in- 
ner accretion disk and the binary orbital plane fjlvanov et al. Ill999l : iLiu 1 12004| ). When a 



SMBHB becomes coalecing due to gravitational wave radiation, the interaction between the 
secondary black hole and the aligned accretion disk would remove the inner disk region and 
l eave a truncated outer accretion disk, leading to the interruption of jet formation in DDRGs 
(ILiu et al. 112003) and to the formation of delay ed X-ray afterglow of a gravitational wave 
radiation burst (IMilosavljevic fc Phinney II2005I ). 

The S- or Z-shaped radio morphological symmetry has been observed f or a high fraction 



of AG Ns and was suggested as an observational evidence of SMBHB by iBegelman et al. 



(|1980[). It is expla i ned with jet precession of periods about between 10^ and 10^ yr (e.g. 



Gower et al. 1982 : Hutchings et al. 19881 : IPunn et al. 2006 ). because of geodetic preces- 



sion of the spin axis of the primary rotating S MBH misaligned with binary total angular 
momentum ( Begelman et al. 1980; Roos lE^SsI ) . the orbital motio n of jet ejecting black hole. 



the d isk precession tidally perturbed by the secondary black hole (IKatz 1119971 : iDunn et al. 



2005 



2006 ), or the precession of a inner warped disk owing to Bardeen- Peterson (e.g. 



Lu fc Zhou 



Caproni et al. II2006I ). In all the models reviewed above, it has implicitly assumed 



that the precession of jet orientation follows the precession of the spin axis of the emitting 
rotating black hole and the rotating axis of the inner region of accretion disk around the 
black hole. Although there are many models for the jet formations, it is usually believed 
that jets form in the inner disk region along the spin axis of black hole or the rotating axis 
of the inner region of accretion disk, depending on the driving energy res ources. Because the 



small characteristic size of jet production region (e.g. iMeier et al. Il200ll ) and the alignment 
of rotating black hole and th e inner region of accretion disk due to Bardeen- Peterson effect 
(IBardeen fc Petterson Ill975l ). it seems reasonable to assume that jet would always orient 
along the rotating axis of both the black hole and the inner regions of accretion disk irre- 
spective of the driving mechanisms. Thus, in this paper we will take the same assumption 
that jets would, if present, precess with the rotating axis of the emitting central black hole 
and the accretion disk. With the assumption, all the models in the literature can explain 
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the observed jet precession, though the results depends on many parameters with very large 
uncertainties. One of the very important question is how to tell which model is the right one 
and to determine parameters. With the improvement of observational instrument, jet pre- 
cession can be observed up to many cycles, which make it possible to measure the precession 
timescale with very high accuracy. In this paper, we investigate the possibility to measure 
the time derivative of jet precession timescale. In the precession models, a rigid-body like 
disk precession is assumed in the literature. Here we take the same assumption. 

A circumbinary accretion disk could be warpe d by a ma s sive SMBHB with r andom 



orbital inclination angle relative to accretion disk (iLiu I |2004J : llvanov et al. Ml999l ). The 
interaction quickly realigns the inner warped disk region and finally the central rotating 
SMBH with the binary orbital plane, while the outer unperturbed disk region far from 
the binary orbit remains coplanar with the galactic plane. This scenario predicates the 
formations of X-shaped radio feature in FRII r adio galax ies and a random distribution of jet 



orientation with respect to the galactic plane (ILiu 1120041 ). A warped disk precesses, leading 



to jet precession. Therefore, before discussing the variations of jet precession timescale, we 
calculate the precession of a warped circumbinary non-massive disk in this paper. Although 
the detailed SMBHB models for jet precession are different, all of them predicates a increase 
of precession tim escale with binary evolution. As SMBHBs in galactic nuclei never get softer 



(iQuinlan Ill996l ) , the secondary black hole always migrates toward the binary mass center and 
the jet precession is expected to be accelerated in the SMBHB models. After calculating the 
jet precession timescale and SMBHB evolution timescale in different models, we show that the 
acceleration of jet precession could reach 20 % or even higher depending on the parameters 
of SMBHB systems and accretion disks. Differential measurement of jet precession can be 
used to distinguish the different precession models and determine the evolution timescale of 
SMBHB in galactic nuclei. With the measurement of precession acceleration, we could also 
determine the kinematic viscosity coefficient of accretion disk and the binary parameters. 

Following the different physical mechanisms driving the evolution of SMBHB in galactic 
nuclei, in Section [2] we start our calculations of the hardening rate of SMBHB and the ac- 
celeration of jet precession with the regime when stellar dynamic friction affects the merger. 
In Section [31 we calculate the evolution of SMBHB and the time variation of jet precession 
because of the interaction of SMBHB with a massive circumbinary accretion disk, which 
is followed by the calculations for the scenario in which the evolution of SMBHB is dom- 
inated by the interaction between SMBHB and a non-massive circumbinary accretion disk 
in Section HI In Section 0, we estimate the acceleration of jet precession because of the 
rapid evolution of SMBHB dominated by gravitational wave radiations. As an example, in 
section [6] we discuss the differential observations of jet precession in a recent merged system, 
NGC1275 (3C84), and the implications to the SMBHB in the object. Our discussions and 
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conclusions on the results are given in Section [71 

2. Hardening of SMBHB because of galactic dynamic friction 

2.1. Unbound SMBHBs 

Two SMBHs in merging galaxies are unbound until the gravitational force between the 
two SMBHs dominates the orbital motion when the separation of SMBHB a is 



GiM + m) . , 

a > an = ^ ^ ~ 1.12 x IOVg 1 + q) 



\ -2 



^200Km/s 



(lYu II2002I ). where o is the one-dimensional velocity dispersions of the primary galaxy, re = 
2GMj(? is the Schwarzschild radius of the primary SMBH of mass M, and q = m/M is the 
mass ratio of the secondary (of mass m) and the primary SMBHs. For a > oh, the evolution 
of SMBHB is dominated by galactic dynamic friction and the evolution timescale of Ta is 
approximately proportional to the separation a 

" rj^\, (2) 



(da/dt) yCLR 

where th = —clr/ Vdy is the dynamic timescale at oh. A negative s ign for the def inition 



of Ta is used because a SMBHB at galactic center never gets softer (iQuinlan Ill996l ). The 



hardening rate of SMBHB due to galactic dynamic friction v^y is approximately 

Vdy ^ -0.l51al/(j^\nA (3) 



Merritt Il2000l ). Here InA ^ 2 is the Column logarithm and Ug is the one-dimensional 



(e.g. 

velocity dispersion of the smaller (secondary) galaxy. Applying the empiric al relation of cen- 



tral b lack hole mass M and the stellar velocity dispersion a of host galaxy (ITremaine et al 



20021) 



log(M/M0) = 8.13 + 4.02 lg(a/200Km/s) (4) 

to both the primary and the secondary galaxies, we obtain th — 8.71x 10^Mg ■'^^'^^■'^^g~^/^■°^(H- 
g) yr and the SMBHB evolution timescale 

ra = -- - 3.73 X lO^Ml-^/^-'glf (^) yr, (5) 

Vdy V 10 / 

where g_i = g/0.1 and M = Mg x IO^Mq. 
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When a > an, jets, if present, would precess because of the orbital motion of the emitting 
black hole with the orbital period 



orb 



27T 



.1/2 



_G{M + m + M,) 



3.93 X lO^Ms 



IOVg 



3/2 



1/2 



yr, (6) 



where M^, is the mass of the stellar objects inside SMBHB orbit and > M + m. Here a 
typical mass M^, ~ 5M is used, because we are interested in a SMBHB with a > oh- This 
is the shortest precession period in a SMBHB system with a given binary separation a and 
has been introduced to explain the helic al jet morphology on pc -scale and periodic optical 



outbursts observed in some blaza rs fe.g. IVillata fc Raiteri Ill999l ) and jet precession at Kpc 



scale or larger in some AGNs (e.g. Wirth et all 19821 ). From Equation ([6]), we can obtain the 



acceleration of jet precession in binary orbital motion because of the hardening of SMBHB 



dR 



orb 



dt 



3 Porh 



(7) 



From equations IQ, ([5]), and ([7]), we get the acceleration of jet precession because of 
binary evolution 



dP 



orb 



dt 



-0.16 



1/2 



1/2 



, .1/4.02 3/4.02 

^yj-s y-i 



10Vg7 \l + q + ^ 



As the precession period is within the observable range of jet precession in the literature, 
equation ([8]) implies that one can measure the evolution of SMBHB owing to dynamic friction 
by detecting the acceleration of jet precession. 



2.2. Bound SMBHBs 



When a < an, a SMBHB becomes bound, while a bound SMBHB becomes hard at a 
separation 



GmM 



Oh 



4o-2(m + M) 
3.260 X IOVgMc 



-1/2.01 Q-l 



l + q 



(9) 



(iQuinlan Ill996l ). The evolution of SMBHB with separation < a < an is still dominated 
by galactic dynamic friction, but the hardening timescale is approximated with r^y oc a'^""'^ 
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( lYu Il2002[ ). where 7 is a fitting parameter in the Nuker law for the inner surface brightness 
profiles of galaxies 



Hr] 



/3-7 / r 
2— Jb - 



r 

1+ I - 



0-1 

V 



(10) 



and 7], (3, Jb, and rb are fitting parameters, too. The break radius rb is the point of maximum 
curvature in log-log coordinates and 7 is the asymptotic logarithmic slope inside rb. For core 
galaxies, 7 < 0.3, while for power law galaxies 7 > 0.5. Therefore, the evolution timescale 
of a bound SMBHB is approximately 



a 

an 



7-0.5 



(11) 



This equation gives 



1.75 X 10^ X 13.03-^M(^+°-°^)/^-°^g:f 
7-0.5 

a-5-7 ■ 



a 



lO^rc 



yr 



(12) 



for Oh < a ^ Ah- Equation (fTTj) implies that the evolution timescale of SMBHB decreases 
with binary separation for core galaxies but slowly increases for power law galaxies. 

If the orbital motion and the angular momentum is dominated by the total mass of 
SMBHB, jet may precess because of the geode tic precession of the sp in axis of the primary 
black hole about the total angular momentum (IBegelman et al. Ill980l ). of the orbital motion 
of the emitting black hole, and of the accr etion disk p recession due to tidal force of the 
inclined secondary SMBH outside the disk (iKatz 1119971 ). If the rotating primary SMBH is 
misaligned with the binary total angular momentum, the spin axis of the primary black hole 
undergoes geodetic precession about the total angular momentum with a period 



P ~ 

geo — 



2.6 X lO^Mggz! 



IOVg 



5/2 



yr 



(13) 



(IBegelman et al. Ill980l : IRoos Ill988l ). This equation gives an acceleration of jet precession 
due to the hardening of SMBHB 



dR 



geo 



5 P. 



geo 



dt 



(14) 



Equations ([I3D, ([H]), and ([I2D shows that 



dP, 



dt 



a \ 



3-7 



(15) 
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for < a < an, implying that the acceleration of jet precession in the geodetic precession 
is too fast to be detected. 

If the orbital plane of a SMBHB is inclined with respect to an accretion disk of radius 
i?d inside the binary orbit, the disk precesses like a rigid body owing to the tidal force of the 
secondary with a precession period 



' '''' + 

lOVcy VlO^^G/ g-icos^ 



-td^5.1xl0VrM|(-£-) (-^) '-^—^ (16) 



( iKatz I I1997I ). where 9 is the tilt angle of the disk plane and the binary orbital angular 
momentum. The disk size could be the total radius extent of an accretion dis k or the 
Bardeen- Peterson radius (IBardeen fc Petterson Ill975l : iNatarajan fc Pringle Ill998l ). When 
the SMBHB becomes hardening, equation (|T6l) suggests that the jet precession will be ac- 
celerated with 

dP^ = _3:^. (17) 

dt Ta 

To obtain equation (|T7j) . we have assumed that the change of the disk radius i?d is insignif- 
icant, comparing to the variation of the binary separation. From equations (|T6l) . ( ITTj) . and 
f[T2|) . the jet precession because of binary-disk tidal interaction is accelerated with 

^^td ^ -0.087 xl3Wf°^-^)/^-°^g:}-°'/^-°'(l + g)^-^ 



dt 

3-5-7 / p \ -3/2 

COS~^ 9. 



The acceleration dPtd/dt of jet precession is very significant for power law galaxies but 
moderate for core galaxies. 

From equations ([6]), ([Tj), and f|T2l) the jet precession period due to binary orbital motion 
will change with time 

^ ^ -7.5 X 10-Mr-)/-^g!r (1 + g)--- , (19) 

where we have taken <C M + m. Equation f[T^ shows that the change may be too small 
to be detectable. 



3. Evolution of SMBHB because of interaction vi^ith massive disk 



When a SMBHB becomes hard at a ~ ah, the evolution timescale Xa may be larger than 
the Hubble time and the binary may stall, if three-body interaction between SMBHB and 
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stellar objects dominates the binary evolution (iQuinlan I Il996l : lYu 1 12002| ) . For a SMBHB 
stalling at the hard radius a^, equations f|T3|l . f|T6|l . and ([6]) together with equation ([9]) give, 
respectively, the constant precession timescale 



P ~ 

geo — 



4.99 X 10«iVV°-^^/^-°^g!/^(l + g)^^/Vr, 



1.77 X lO^M 



1.02/2.01 (1-1 



i?d \-"' 1 



l + g)5/2 \^io4rGy' cos^'^^' 



3/2 



^ 1.63xl0X-°'/'-°'7T^yr. 

(1 + g)^ 



(20) 
(21) 

(22) 



However, gas disk exists at the central region of AGNs, which should interact with 
SMBHB. In the AGN unification model, the size of accretion disk around central SMBH is 
order of lO^rc while the broad emission line region and thick dust torus outside the accretion 
disk can be as large as ~ lO^rc ~ 10 pc. Because the dust torus is geometrically thick 
and massive, the interaction between the secondary black hole and the dust torus is linear 
and the secondary cannot open a gap, probably leading to a rapid ty pe I migration of the 



secondary toward the mass center (e.g. iPapaloizou &: Terquem 1 12006| ) . If the accretion disk 



is geometrically thin and coplanar with the binary orbital plane, the secondary SMBH with 
mass ratio q > gmin = ^aS"^ ~ 3.2 x 10~^a_i (5/0.01)^ will open a gap in the accretion disk 
and exchanges a ngular momentum with disk gas via non-linear Lindbl ad resonant binary- 
disk interaction (ILin fc Papaloizou Ill986l : lArmitage fc Natarajan 1120021 ). Here, the viscous 
parameter a = 0.1a_i is defined with the shear viscosity in r-0 plane, i>i = aCgH with H 
the scale height of the unperturbed accretion disk and Cg the sound speed, and S = H/r 
is the disk openning angle at radius r. The migration of the secondary black hole is called 
Type II migration. If the disk mass inside the binary orbit is larger than the mass of 
the second ary SMBH, the migration timescale of the secondary SMBH is the disk viscous 



times cale (ILin &: Papaloizou 1 1 19861 : lArmitage &: Natarajan I l2002l : IPapaloizou &: Terquem 
20061 ). If the circumbinary disk is massive, the secondary SMBH migrate also on a disk 



visco us timescale even if the orbital plane and the accretion disk is misaligned (llvanov et al 



19981 ). Because the mass rat i o of a SM BHB formed in galaxy mergers within Hubble time 



is g > 10 ^ > gmin (IYu II2OO2I : iLiu II2004I ). we consider only type II migration. 



For a type II migration, the secondary black hole migrates inwards on a viscous timescale 

(23) 

From equations 0231) and Ui = aCgH, we have 



a 2a2 

^ -f- 

'a ~ ''i/ — — r, 

f r 3 Ui 



Ta - 1.66 X lO^MgaZldzl 



5/4 



yr 



(24) 
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where a_i = a/0.1 and S-2 = So/0.01. Here for convenience, we have written 



'^''[w^o) • '''' 

where is the disk opening angle at r = IO^tq and weakly depends on a, the central 
black hole mass M, and the accretion rate m. For a standard a-disk, A = 1/8 if the 
disk is gas pressure and free- free absorption dominated, while A = 1/20 if the disk is gas 



pressure and electron scattering dominated (jKato et al. Ill998l ). To get equation ( IMl) . we 



have assumed that the disk is gas pressure and free-free absorption dominated with A = 1/8, 
which would be valid for r > 2.6 x lO^rQmil faato et al. 1998 ). m = M/MEdd = 0.1rh_i 



is the dimensionless accretion rate, and the Eddington accretion rate Msdd = -Z^Edd/O.lc^ is 
related to the Eddington luminosity LEdd = 1-26 x lO^^Mgerg s~^. 

Because the disk is massive, the total angular momentum of the binary-disk system is 
dominated by the disk mass and there is no geodetic jet precession around binary orbital 
angular momentum. However, the inner disk region misaligned with a rotating central SMBH 
may be warped and become aligned due to Bardeen-Peterson effect ( iBardeen &: Petterson 



19751 ). jets may precess because of the tidal interaction of a misaligned secondary black hole 



and the warped inner disk. From equation (|T6l) . ( JT71) . and (IMl) . we have 

^^^^ - -2.QlM,a.,5\ ( ( (1 + qf' qZl cos'' 6, (26) 



dt • o 1 VIOVg/ \50r, 

where tbp with typical value of order ~ SOrc is the Bardeen- Peterson radius out to which 



the accretion disk flow is a ligned with the black hole spin axis ( iBardeen fc Petterson Ill975 



Natarajan fc Pringle Ill998l ). Equation fl26l) implies that the acceleration of the jet precession 



is significant and can be detected very easily. 

Jets will also precess because of the orbital motion of the emitting primary black hole in 
the case of massive circumbinary accretion disk. However, the acceleration of jet precession 
due to binary orbital motion 



P / ^ \ 1/4 

"-'^orb <J ^orb 



. , - 2.5X10-W12 Ut^^ (l + q)-'^' (27) 

at 2 Ta VIO^G/ 

may be too small to be detectable. 

When the inner disk region becomes aligned with the rotating black hole but misaligned 
with the outer inclined accretion disk, the misaligned disk region and the spin axis of central 
rotating black hole would precess with a precession timescale 



BP 



1.51 X 10'al'/V('mZ'/'M~'/'' yr (28) 
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( iNatarajan fc Pringle Ill998f ). where a* is the spin parameter of the primary SMBH. Equa- 
tion fl28|) shows that the jet precession due to Bardeen-Peterson effect is independent of the 
evolution of SMBHB and does not change with time on the timescale that we are interested. 



Evolution of SMBHB because of interaction with a non-massive disk 



When the secondary black hole migrate inwards on a viscous timescale and reaches a 
critical radius r^, the disk mass inside the binary orbit will equal to the mass of the secondary 
black hole. In a gas-pressure and electron-scattering dominated a-disk, the unperturbed disk 
surface density is 



2.4 X lO^a' 



4/5^1/5^3/5^-3/5 



1 



m 1 r 



1 ' 3 



g cm 



(29) 

( iKato et al. Ill998l ). where rs = t/I^^tq. Note that we have used different a prescription. 
From equation fl29l) . we can estimate the disk mass inside radius r 



Md ~ yvrSr^ ~ 4.9 x IQ^aZl'^M^ 
When Md = m, from equation fl30|) we have 

= 8.6 X lQ\''lla1lm-_TM,"\G. 



11/5 . 3/5 7/5 

m \ r 



1 '3 



(30) 



(31) 



When a < Tm, the disk mass Md inside the binary orbit is smal ler than the mass of 



the secondary and the migration of the secondary will be reduced ( ISyer &: Clarke I Il995 



Ivanov et al. Ill999l ). If the disk is gas-pressure and electron-scattering dominated, the mi- 
gration timescale is approximately 

X5 \ 5/19 16/19 / ;,,r.\ -14/19 

152 J 



152 



112 



16 
IT 



Md 
m 



(32) 



( llvanov et al. Ill999l ). where Md = Mt^y is disk mass inside the binary orbit and t^, is the 
viscous timescale of unperturbed accretion disk at r = a. Taking ti, = (2/3)a^/z/i and from 
equations (l25l) and (l32l) . we have 

7/19 



8.95 X 



14/19^14/19j^5/19^ -5/19^^-10/19 



1 



Q-i 



a 



-1 



'5 



-2 



lO^rc 



yr 



(33) 



for a gas-pressure and electron-scattering dominated disk with A = 1/20. 



If the rotating primary black hole is inclined to the binary total angular momentum, 
its spin axis will precess geodeticaly. From equations (fT3!) . (fT^ . and (!33|) . we obtain the 
acceleration of geodetic precession 

81/38 

(34) 



geo 



dt 



-2.3 X 10-WTM's'^''qZf''a'/l'5'_T 



IOVg 
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When the disk mass inside the binary orbit is less than the mass of the secondary 
black hole , the secondary will warp, twist, and realign the inner accretion disk on a short 
timescale ( llvanov et al. Ill999l ). However, the realignment will stop at the Bardeen- Peterson 
radius tbp and the disk region at r < tbp will rem ain aligned with rotating primary black 
hole and misaligned with the binary orbital plane (ILiu II2004I ). So, the inner accretion disk 
has r < tbp and thus the jet orientation precess due to the tidal interaction of the secondary 
black hole. Equations (fT6ll . (fT7|l . and (!33|) give the acceleration of the precession period 
because of binary-disk tidal interaction 



td 



-4.8 X lQ-'rn'TMf''\-_f''\T5'T 



. 50/19 , 

a \ ' /rep 



- — i — i ^103^^ 

The time variation of jet precession is very small. 



SOrc 



'3/2 



} + 



1/2 



COS 6' 



(35) 



When the disk mass Md within the binary orbit is less than that of the secondary black 
hole, a SMBHB with an inclined orbital plane war ps and realigns the inner disk region 



outside its orbit to a typical transitional radius ( llvanov et al. Ill999l ). The transitional 



radius rai of the inner warped and the outer unperturbed disk regions depends on how warps 
communicate in the di s k. Ta king into account the internal hydrodynamics of the disk itself, 
Papaloizou &: Pringle I (119831 ) showed that for a > 6 = H/r warp transfers on a timescale 
t^p ~ 2r^/3^'2, where z/2 is the vertical viscosity and relat es to the shear viscosity z /i in r-(j ) 
plane with i>2 ~ i'] f n/ (2q;^) and /o = (1 + 7a^)/(l + a^/4) (IKumar &: Pringle Ill985l : Kumar 
199dlOgilvie]ll999[ ). 



The quadrupole contribution of the secondary black hole to the gravitational potential 
would causes the precession of the major axis of an elliptical orbit in the disk with frequency 

3 



n 



ap 



(llvanov et al. 1 119991 ) . where Qk is the Kepl e rian a ngular velocity at 



precesses with frequency Q^p = — ^ap- iLiu I (120041 ) showed that for /( 
radius rai can be estimated by using t 



wp 



fla_p, which gives 



(36) 

The lines of nodes 
= 1 the transitional 



(37) 



where r^i < r^. The precession period Pap of the aligned inner disk is determined by the 
precession of the lines of nodes at Tai 



P. 



ap 




^q21A/2(1+A) 

V C 



a 
01 



7/4(l+A) 
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(3-4A)/2{l+A) 

^(3-4A)/4(l+A) (38) 

For a gas-pressure and electron-scattering dominated thin disk, A = 1/20 and the precession 
period is 

P.P ^ 2.52 X 10^ yr Msq%\%' f,~'^'6Zl'^' ( ) . (39) 



IOVg/ 

Because the warp transfer timescale t^p = (2a//o)^v is much smaller than the viscous 
timescale for a standard thin disk with a ^ 1, an assumption of rigid-body like disk 
precession is reasonable. If the primary SMBH is rotating and misaligned with the binary 
orbit plane, warp transfers quickly inwards and stall at the Bardeen-Peterson radius rep 



(IBardeen fc Petterson Ill975l ). The disk region within rep and the spin axis of the rotating 



primary BH would also precess with a timescale Pbp- 

From equation fl38l) . the acceleration of the precession of a warped circumbinary disk is 

d In Pao dlnM 7 dln6o 3 — 4A c? In a , , 

— = H (40) 

dt dt 2(1 + A) dt 2(1 + A) dt ^ ' 

where we have assumed a constant mass ratio q during the evolution of SMBHB. In accretion 
disk theory, the opening angle 5q depends on accretion rate and the mass of central black 
hole 

(5o oc m^'M^ (41) 



with /i > and C > (IKato et al. Ill998l ). Substituting equation fHTl) into equation (j40|) . we 
have 

rIP 7// P R-A\ P 

(42) 



C^Pap 7/i Pap 3-4A Pap 



dt 2(1 + A) 2(1 + A) Ta ' 

where = —rn/{dfn/dt) is the variation time scale of accretion rate, which is equivalent to 
the typical lifetime of the AGN and may be determined by the environment, for example, 
the supply of the gas from the galactic disk to the accretion disk and the interaction of 
accretion disk and the stellar objects passing through the disk. To obtain equation (142|) . 
we have assumed that the mass of the primary SMBH is insignificant on the timescale 
which we are interested in here, implying that the variation timescale of the black hole mass 
tm = M/{dM/dt) is much longer than the binary hardening timescale t^. Equation (H2l) 
suggests that the decrease of accretion rate decelerates the jet precession but the hardening 
of SMBHB accelerates it. If a SMBHB is long-lived and passes through the active phase of 
a galaxy, namely <C Ta, equation (l42l) gives a deceleration rate of jet precession 



dt 3 
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for a gas-pressure and electron-scattering dominated standard thin disk with /i = 1/5 and 
A = 1/20. Equation fH3l) suggests that from the measurement of jet precession and its 
deceleration rate, we can determine the disk evolution and the life time of an individual 
radio source. For a short-lived SMBHB with ^ r^, jet precession is accelerated 

for a gas-pressure and electron-scattering dominated thin disk. 

From equations (jSl) and ( l33l) . we obtain the acceleration of the precession of a warped 
disk due to the reduced migration of the secondary black hole 

^ ^ -3.8 X lO~WTMl'/''qZt^''fo'^'alT'SZl''^'' {l^) ' ^^^^ 

/ \ 3/2 

At a ~ IOVg, the orbital period of a SMBHB is Porb ^ 8.8 f yp^J {i+q^n and the 

jet precession because of the binary orbital motion would be nearly constant as dPorh/dt ^ 
3 fprb < 1 n-5 

2 Ta ~ ^ ■ 



5. Rapid Evolution of SMBHB owing to gravitational wave radiation 



When a is order of lO^rc, the loss of the orbital angular mornentum because of gravita- 
tional wave radiation becomes important ( lArmitage fc Natarajan 1120021 ) and the in-spiraling 
velocity of the secondary black hole due to gravitational wave radiation is 

UG^M^q (1 + q) 



"gw 



-4.8 X 10" 



lO^G 



/g_i (1 -|- q) cm s~^ 



( jPeters fc Mathews Ill963l ). where / is a function of eccentricity e 

(1-eT^^ 



/ 



, 73 2 37 4 

1 H H 

24 96 



(46) 



(47) 



At large a, the evolutions of the SMBHB and the accretion disk are coupled and the migration 
timescale of the secondary SMBH is given by equation ([321) • When a is small, the loss of the 
binary orbital angular momentum is dominated by the gravitational wave radiation and the 
hardening timescale due to the gravitational wave radiation is given with 

ra = -a/agw^l.95xl0^M8g:J(l + g)-V"'(Y0^) ^^^^ 
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At a critical radius a = Ogw, the in-spiraling timescale due to the gravitational wave radiation, 
Tgw, approximately equals to the migration timescale of the secondary SMBH because of the 
interaction between the SMBHB and a non-massive accretion disk, Ta. From equation ([3SD 
and = — flgw/flgw, we obtain 

^ 4.27 X 10'rom-_l'/''M-''/''q'}l'' (1 + g)^^/^^ aZf'SzT'r^''. (49) 



From equations (fT3!) . (IT^ . and (l48ll . we obtain the acceleration of jet precession because 
of the geodetic precession 

%.3.3xl0-(l.,)/(^)"^ (60) 

Equation ( l50l) implies that the time variation of jet geodetic precession depends only on the 
binary separation and is independent of the parameters of accretion disk and the SMBHB. 

From equations IQ, ([7]), f|T6l) . f|T7|) . and fHHj) . the accelerations of the jet precession 
owing to the tidal interaction of binary-disk and to the orbital motion are, respectively. 



td 



dt 



(51) 



^ . 2.1xlO-V.(l + #/V(^)'"^ (52) 
The acceleration are insignificant. 



From equations fl39|) . fl44j) . and fl48l) . we get the acceleration of jet precession because of 
precession of a warped circumbinary disk 

^ ^ -0.80g!/^(l + q)fa%'6Zl'^'f-'/' . (53) 



6. Differential jet precession and SMBHB in NGC1275 

In previous sections, we discussed the evolution of SMBHB in different driving regimes 
and the corresponding acceleration of jet precession. Equations ( 1141) . (ITTll . (JTj), and (14^ 
suggest that the acceleration of jet precession because of SMBHB hardening can be written 
integrally 

dP P 

^ ~ -A^, (54) 
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with I < A < 3. Our calculation show that the acceleration of jet precession because of 
SMBHB evolution is significant and could be measured on the timescale of jet precession 
period. If we measure the jet precession period Ppr and compute the acceleration rate 
dPp^/dt, we can determine the evolution timescale of a SMBHB in galactic nuclei with 

ra ^ -A^. (55) 

pr 

With these calculations, we will discuss, as an example, the differential jet precession in the 
radio galaxy NGC1275. 



6.1. Jet precession with constant timescale 



The FRI radio galaxy N GC1275 (3C84) is a recent merger system at redshift z = 0.01756 



[e.g. iHoltzman et al. I Il992l ) and the mass of central SMBH is measured w ith molecular 
gas hydrodynamic method to be M = 3.4 x 10^ ( iWilman et al. II2005I ). The object 
has a bol ometric luminosity Lhni — 1-07 x 10^^ ergs s~^ (with Hq = 75kms~^ Mpc~^ and 



go = 0.5) (IMarchesini et al. \\2004 ). With the measured black hole mass and the bolometric 



luminosity, we obtain the dimensionless accretion rate m ~ 2.4x 10~'^elp iDunn et al. I ( 20061 ) 
imaged the S-symmetrical morphologies of jets and emission line structure and differentially 
measured jet precession timescale by identifying four components of different orientations 
in order of formation: ancient bubbles, ghost bubbles, outer lobes, and inner jets. The 
observations of precession angle A0, time difference At, the number of complete cycles n 
between two comp onents in succession, the observed precession period Ppr are taken from 
Dunn et al. I ( 120061 ) and summarized in Tabled! The number n is the precession cycle number 
between two successive components and is estimated with Ppr = 36O°At/(A0 + n360°). 
The observations suggests that the activity of NGC1275 is intermittent and the precession 
timescale are significantly different for different episodic activities. Intermittence of activity 
and significant differences of jet precession periods for differ eii t activity episodes are also 



observed in the Seyfert 1.5 galaxy Mrk 6 (iKharb et al. II2006I ). For both NGC1275 and 
Mrk6, the observations shows that the jet orientation at the beginning of each episodic 
activity are significantly different from that at the end of last episodic activity, implying 
that the spin axis of central black hole precesses even when the source is dormant or at very 
weak activity. 



Dunn et al. I (120061 ) assumed that the jet precession in NGC1275 remains steady for 



many cycles and the observed differences of precession timescale are due to the missing of 
different precession cycles when no bubble detaches. With the assumption, the significantly 
different precession timescale from Ancient bubbles through Ghost bubbles to outer lobes 
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Table 1: Observations of jet precession in NGC1275 (3C84) bv lOunn et al. I (120061 ) . Both 
At and Ppj. are in units of 10''yr. 

North South 





A0 




61° 


A0 = 


- 113° 

— X X fj 








3 8 


At - 


-72 


n = 


P 

^ pr 
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P = 

^ pr 


= 22.9 


n = 1,2 


P 

pr 




3.25 


P - 
^ pr 


= 3.11 
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A0 
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A0 = 


= 293° 




At 




6.1 


At = 


= 6.2 


n = 


p 

^ pr 




8.01 


p = 

^ pr 


= 7.62 


n = 1 


P 

pr 




3.46 


P = 
^ pr 


= 3.42 


Outer— >Jet 


A0 




201° 


A0 = 


= 249° 




At 
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At = 


= 1.9 


n = 


p 

pr 
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p - 

pr 


= 2.75 


Ancient— >Outer 












n = 


{dPpr/dt) 




-2.55 








(Ppr) 




15.23 






n = l,2 






3.31 
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n=0 


{dPpr/dt) 




-1.30 












5.27 






n=l 


(dPpr/dt) 




-0.19 








(Ppr) 




3.08 







- 18 - 



are reconciled with one period Ppr — (3.31 ± 0.46) x 10^ yr, but the curr ently active jets sti ll 



precess with a significantly shorter timescale Ppr — (2.72 ±0.54) x 10 yr (iDunn et al. II2006I ). 



Here the errors have included the measurement error of precession timescale, cp/Ppr ~ 20%. 
The cycle number is given in Table [H The jet precession from Ghost bubbles to active 
jets is accelerated with (dPpr/dt) ~ —0.19, which implies that the jet precession is due to 
SMBHB at center. Equation (!55|) gives a model-independent evolution timescale of SMBHB 
in NGC1275 r^h — 3.62 x 10^ yr. If we knows the binary separation a, we can calculate 
a model-independent binary hardening rate or the migration velocity of the secondary Vob — 
-54 (ol^) (f cm/s with 4/3 < A < 3. 

Jets precess with constant timescale through several duty cycles of activity implies 
that the precession in NGC1275 is independent of the accretion. All the models for jet 
precession depending on accretion disk are excluded and the only reasonable scenarios for 
the jet precession from one bubble to another are the geodetic precession or binary orbital 
motion. If the precession is due to the orbital motion of a SMBHB, the observed period 
and equation ([6]) gives Oorb — 1-81 x IOVq [(l + g + ^) /^Y^^ ^ 5.9 x lO^pc, which is 
much larger than the bound radius oh — 7.1 x lO^rn and implies an unbound SMBHB 
in NGC1275. An unbound SMBHB is consistent with the observations of recent merger. 
However, equation ([H]) suggests a variation of precession timescale from the ancient bubbles 
to outer lobes for the south components 



APpr (iPorb 
Ppr ~ dt 



2P 

pr 



1 + ? + ^ 



which together with the assumption of steady precession gives an upper limit q ^ 2 x 
10^^ [(l + g + ^) /5] Equation([5]) shows that to form such a binary in a minor merger, 

— iqT^ ) yi') where M^, 

is the stellar mass within the binary orbit with a ~ 20Kpc. A minor merger with q 10^^ 
is unlikely to be observable for such a long timescale. In this scenario, a steady precession 
from the ancient bubbles to the outer lobes is also inconsistent with the acceleration of jet 
precession from the Ghost bubbles to the present active jets. 

The second possible precession scenario independent of accretion rate is the geodetic 
precession of the primary SMBH. From equation ( |T3|) . the constant precession timescale 
implies 

ageo - 6.7 X Wr^q^ll (57) 

which is smaller than Oh for q > 2.0 x 10^^. When the source is at very weak activity or 
dormant with an accretion rate much smaller than its current accretion rate and also the typ- 
ical accretion rate for FRI radio galaxies, namely rh <^ 10"'^, the accretion disk cannot be a 
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standard thin di sk but geometrically thick and optically thin advection dominated accretion 



flows (ADAFs) (INarayan &: Yi Ill994l : iMeyer fc Meyer-Hofmeister Ill994l : lAbramowicz et al 



19951 ). The interaction of an AD AF and the secondary black hole in a binary system is 
dynamically negligible (iLiu 1 12004| ). For such a binary-disk system, the total angular mo- 
mentum is dominated by the orbital angular momentum and the jet precession is geodetic. 
For a ADAF-binary system, the jet precession other than the geodetic precession is the 
orbital motion with period Porb — 5 x 10^ yr. For an accretion disk with accretion rate 
fn = 2.4 X 10~^, the inner region is ADAF and the outer part of the disk is a standard thin 
disk. The transition radius between the two different accretion modes is 



3.1 X 10>G 



(5^ 



(Meyer et al. 


2000: 


Liu et al. 


2002) 



1.4 X 10 ^ but less than the transitional radius r 



2x 10-^a 



3 -20/11 ■ 15/11 



m 



Equations fl57j) and fl58l) shows a^co > rtr for q > 

o 1 in4 5/7 4/7 . -3/7 

2.1 X 10 rGg_i a_i m ' 

3 



-1 



for q 



> 



where m_3 = m/10 . The secondary SMBH migrates inwards owing 
to the interaction with the standard accretion disk on a timescale given with equation fl55]) 



Ta ~ 3.88 X lO'gr/ 



8^84/95^-5/19^-10/19 



a 



(59) 



From equation fl34p . the time derivative of geodetic precession because of the migration of 
the secondary interacting with a circumbinary accretion disk is 



rfP 



geo 



dt 



-0.21g_° 



-84/95 5/19 el0/19 



a 



-1 



-2 



(60) 



The observed acceleration {dP^r/dt) ~ —0.19 and equation fl60|) suggest a binary mass ratio 
q ^ 0.11ci;^Y^^'^-2''^^ the secondary has mass m ^ 3.8 x 1O''M0q;^Y^'^'^-2'''^^- The binary 
separation is about a ^ 7.0 x lO^oa^''^'^^^^^^"'^ ~ 0.23 pc and a < Oh ^ 1.8 x lO'^rca'^^i^^S'^'z^'^ ■ 
The results are insensitive to the disk parameters. 



However, the secondary black hole should warp the standard thin disk and the warped 
disk would precess, probably leading to the precession of jet orientation with a timescale 
given with equations (155]) and 



P 

ap — 



1.08xm'q'_{'a'!^f-'/'6Zf'jT 



(61) 



If the precession timescale from the outer lobes to active jets is due to the precession of 
the warped disk, equation fl6Tl) and the measured period Ppr = 2.72 x 10^ yr give q = 
0.21a_'^^^^^ fQ^^^^6'^2^ . The secondary has mass m ~ 7.2 x 10'^ MQa_l^^^^ Jq^^^^ S'^J^ and the 



^1 ^0 "-2 

binary separation is about a 



5/9 f5/9 rlO/9 



25/18 ^25/18^25/9 



9.1 X WrGaZTfo' 



-2 



'©"-1 Jo "-2 

0.29 pc and a < Oh 



3.1 X 







The secondary migrates inwards because of binary-disk interaction. 
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leading to a time variation of jet precession timescale 



^ ^ -1.4xlO-2a^f/o-^■«^r^^ (62) 

dPgco ^,-1 85/57 y.-70/57 r-110/57 



dt 



-Qlla'^T ' ' ■ (63) 



Equations flB^ and (jHSD implies that the time variations of jet precession cannot be detected 
because of the low observational accuracy. Therefore, the different precession timescale be- 
tween the outer lobes and the active jets is most probably because of the different mechanism 
for jet precession. 



6.2. Rapid acceleration of jet precession? 



The argument for a steady precession in the object given by iDunn et al. I (120061 ) is 
that if the precession is speeding up, the acceleration would be very rapid and over the 
courses of about 1.5 rotation the precession timescale changes by around a factor of 10. 
However, our theoretical calculations suggest that a rapid acceleration is possible and there 
is no a priori requirement for constant precession timescale. In this section, we discuss the 
implications of a rapid acceleration of jet precession. We compute the precession timescale 
and the time derivatives for n = in Table [H From the averaged precession timescale from 
the Ancient to the outer lobes (-Ppr) — 15.23 x 10^ yr and the averaged time derivative of 
the precession timescale from the Ancient bubbles to the outer lobes {dP^r/dt) ~ —2.55, 
we have the model-independent evolution timescale of SMBHB — (j^) 8.96 x 10^ yr, 
which is about three times smaller than the timescale obtained with the assumption of a 
steady jet precession from the Ancient bubbles through the Ghost bubble to the outer lobes. 
While from the averaged precession timescale from the Ghost bubbles to the active jets 
(Ppr) ~ 5.27 X 10'' yr and the averaged acceleration of the precession from the Ghost bubbles 
to the active jets {dPpr/dt) ~ —1.30, we compute the model-independent evolution timescale 
of SMBHB Tgj ~ (^) 6.08 x 10^ yr. 

As the precession from the Ancient to the active jets is continuous without interrup- 
tion when the activity of the object varies significantly, the possible mechanisms for jet 
precession are the binary orbital motion and the geodetic precession of the primary. If 
the precession is due to binary orbital motion, the averaged precession timescale from the 
Ancient to the outer lobes (-Ppr) — 15.23 x 10^ yr and equation ([6]) give an averaged bi- 

nary separation Oao — 5.01 x 10 re I — ) — 1.46 Kpc. From the averaged time 



derivative of the precession timescale from the Ancient bubbles through the Ghost bub- 
bles to the outer lobes {dP^-^/dt) ^ —2.55 and equation ([8]), we obtain the binary mass ratio 
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/ ^^ M, N 4.02/9 

Qao ~ 0.92 ( — j . Meanwhile, the averaged precession timescale from the Ghost 
bubbles to the active jets (-Ppr) — 5.27 x 10^ yr and equation ([6]) give an averaged binary 

separation agj ~ 2.47 x 10''rG ( — ) — 0.80 Kpc. From the averaged time derivative 
of the precession timescale from the Ghost bubbles through the outer lobes to the active jets 



1^^x4.02/9 



{dPpr/dt) ~ —1.30 and equation ([8]), we have binary mass ratio ggj ~ 0.60 , ^ 
The two averaged time derivatives give a consistent mass ratio and suggest a major merger 



with (g) ^ 0.76 y — 5"^j • From galactic dynamics (iBinney &: Tremaine 1119871 ). it is 
expected that the evolution timescale of a SMBHB due to dynamic friction is nearly propor- 
tional to the separation, Xa oc a. Our results give Tg^/r^ — 1-5 and aao/ctgj — 2.0, which are 
consistent with the predications very well and give an averaged dynamic friction velocity 

.\ .^4x 106 ( i±i±JVL^ 

dt 



(^) ~ _1.54 X 10^ ( ^^"i^^M I ^^/g (g4) 



where is the stellar mass inside the binary orbit at separation a ~ 1 Kpc. 

The alternative for jet precession independent of source activity and the accretion is the 
binary geodetic precession. From equation (ITSjl . to obtain the averaged precession timescale 
from the Ancient to the outer lobes, we have the binary separation Ogeo — 1-24 x lO^rcQ'^'^i^. 
Because a^^o ^ cth — 1-77 x 10^rGq'_i/(l + g), the binary is hard and the secondary may 
migrate inward when it interacts with a light standard disk, leading to the acceleration of 
jet precession 

^ ^ -QAlQm^TcT'o^'ll'^'T. (65) 

where we have used accretion rate rh ~ 10^'^ and equation (!34l) . Equation (165|) and the 
measured time derivative give g ~ 1.3 x 10^^m^''3 a^Y*^^^^ 2''^^- Similarly, from the observation 
of jet precession from the Ghost bubbles to active jets, we have Ogeo ~ 8.13 x lO^rcg^'^i^ and 
q ~ 0.97 X 10~^mj^3 a^Y^^'^-2''^^- Although the estimated mass ratios are consistent with each 
other within the uncertainties of accretion rate, the hardening of a hard SMBHB depends 
on the accretion and the migration of the secondary should stop when the source becomes 
dormant and the accretion disk becomes ADAF. Even if the migration could happen when 
the source is luminous and forms the bubbles, the time scale to form a hard binary with 
mass ratio q ~ 10^^ is r > 2 x 10^ yr and is inconsistent with the scenario of recent merger. 
Therefore, the scenario of geodetic precession for jet precession is less favorable. 
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7. Discussions and conclusions 

SMBHBs are expected by the hierarchical galaxy formation model and may have been 
observed in many AGNs. Jet precession observed in many AGNs is one of the observational 
evidences. In this paper, we start our work with the discussion of different mechanisms 
for the jet precession, including (1) the geodetic precession of spin axis of central primary 
SMBH around total angular momentum, (2) the orbital motion of the SMBH ejecting plasma 
jets, (3) the inner disk precession because of the tidal interaction of an inclined secondary 
black hole, (4) the precession of a circumbinary disk warped by the SMBHB, and (5) the 
disk precession because of Bardeen- Peterson effect. The precession of a circumbinary disk 
warped by a SMB HB is discusse d first time. We did not discuss the precession model due 



to disc instability (iPringle 1 119971 ) . because it suggests a stochastic precession rather than a 
regular precession and is inconsistent with the observations of jet precession in most AGNs. 
Although Bardeen-Peterson effect does not directly connect to the presence of SMBHB, the 
origin of misalignment between the rotating central black hole and the accretion disk may be 
due to the interaction of accretion disk and an inclined SMBHB. When the inner disk region 
becomes misaligned with the binary orbital plane owing to the Bardeen-Peterson effect, the 
tidal interaction of the secondary to the warped inner disk also leads to jet precession. 

In these scenarios for jet precession, the precession timescale ranges from order of years 
to much longer than 10^ yr, depending on the parameters of SMBHB and accretion disk. 
However, the parameters are very difficult to determine and the observations of jet precession 
timescale cannot give restrict constraints on the models and the parameters, as they are 
degenerate. Therefore, we suggested to observe one more quantity, the time variations of 
jet precession timescale, to resolve the parameters. We calculated the time variation of 
jet precession in different models and showed that jet precession is always accelerated in 
an evolving SMBHB system. The acceleration of jet precession is related to the evolution 
timescale of SMBHB with ^ ~ -A^, resulting from the fact that all SMBHB models 
for jet precession predicate a relation Ppr oc a'^ with A > and that a SMBHB in galactic 
nuclei never gets softer. The parameter A slightly depends on model with 4/3 < A < 3. Our 
investigations also show that jet precession because of Bardeen-Peterson effect is decelerated 
with AGN evolution. Our results suggest that the sign of the time derivative of precession 
timescale can be used to identify SMBHB models from the others. 

Our calculations show that the time variation of jet precession is proportional to the 
timescale ratio of jet precession and SMBHB evolution. We analytically estimated the evo- 
lution timescale of SMBHBs at different evolution stage and the time variation for jet pre- 
cession in different models, based on our current knowledge of SMBHBs. Our calculations 
show that for an un-bound SMBHB the mechanism for jet precession is the orbital motion 
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and the quick binary evolution because of galactic dynamic friction leads to around 20 % or 
higher acceleration rate of jet precession timescale. For a bound SMBHB system, jet preces- 
sion could be due to geodetic precession of the rotating primary black hole, disk precession 
because of tidal interaction between a standard accretion disk and the secondary, and the 
binary orbital motion. At this stage, the evolution timescale of SMBHB depends on the inner 
surface brightness profiles of galaxies and is estimated with an a symptotic analytic relation 



of the binary hardening timescale and the separation given by lYu I (120021 ) . Although the 
estimate is very rough, our results suggest that the evolution timescale of SMBHB is several 
order of magnitude shorter than the geodetic precession timescale and longer than the binary 
orbital period. So, the geodetic precession is not significant and the time variation of orbital 
motion is difficult to measure. However, if the jet precession is due to the tidal interaction 
of the secondary black hole and an inner misaligned accretion disk, the acceleration rate of 
jet precession could be a order of 10 % or higher. 

When a SMBHB becomes hard and stalls, jet precession may be steady for a timescale 
longer than the Hubble time. Becaus e the migration of the se condary SMBH due to the 



interaction with an ADAF is negligible (INarayan ll2000l : lLiu II2004I ). a nearly steady precession 



jet may be possible if a SMBHB interacts with an ADAF or the precession is due to Bardeen- 
Peterson effect. The fundamental difference between the two scenarios is that the accretion 
disk is a geometrically thin standard disk in the later but geometrically thick ADAF in the 
former. The accretion mode depends on the relative accretion rate m, which one could infer 
by estimating the bolometric luminosity and central black hole mass. However, it is most 
probable that a SMBHB interacts with a standard disk either massive or light. We compute 
the time variation of jet precession because of SMBHB-accretion disk interaction and show 
that in both cases the binary-disk interaction would lead to a significant acceleration of jet 
precession: the acceleration is significant for jet precession because of tidal interaction of 
the secondary and a massive disk but both of geodetic precession and warped circumbinary 
disk precession in the case of non-massive disk, depending on the parameters of the binary 
system and the accretion disk. 

When the evolution of a SMBHB is dominated by the gravitational wave radiation, the 
binary separation is about hundreds of Schwarzschild radius or less. If a jet ejects from 
the central black hole, it precesses because of the black hole geodetic precession, of the 
tidal interaction of binary and inner misaligned disk, of binary orbital motion, and of the 
precession of warped circumbinary disk. Our calculations show that the precession of a 
warped circumbinary disk is strongly accelerated owing to the migration of the secondary 
black hole because of gravitational wave radiation. The acceleration of jet precession because 
of geodetic precession is also very significant for a SMBHB with non-zero eccentricity. 
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When we calculate the jet precession timescale and its time variation, we have assumed 
that disk precession is rigid-like and the jet precession is directly related to it. Although 
almost all the disk precession models for jet precession in the literature adopted the same 
assumption and have successfully explained the jet precession in some AGNs and micro- 
quasars (e.g. SS433), this assumption need more discussions. Whether the assumption of 
a rigid body like precession is valid or not depends on the warp transfer in the disk. As 
we have discussed in Section HI the transportation of warps in disk depends on the vertical 
shear viscosity and the transfer timescale at the transition radius between the warped and 
unperturbed disk regions is on the same order of the precession ti mescale both for the 



Bardeen- Peterso n effect and the warped circumbinary light disk (e.g. iNatarajan fc Pringle 



19981 : iLiu II2004I ). Therefore, the rigid body like approximation for disk precession is correct 
on the zero order of magnitude. As the jet precession and its time derivative depend on 
disk characters in a similar way, the relationship of the ratio of the precession timescale and 
its variation rate, Ppr/Ppr = — -Ppr/ (^^^^ {'t^/o)-, and the SMBHB evolution timescale t^, 
would be expected to be insensitive to how warps transfer in the disk. 

Following our theoretical investigation on the acceleration of jet precession, we discussed 
the implications of the differential observations of jet precession in NGC1275 (3C84), a 
recent-merger radio galaxy. The differential jet precession have been measured between four 
different components in order of formation: ancient bubbles, ghost bubbles, outer lobes, and 
the active jets. Between the formation of different components, the activity of the object 
becomes very weak or the source is dormant. The precession ti r aescal e are significantly 



decreased with time among the different components. iDunn et al. I (120061 ) assumed a steady 
jet precession and the acceleration of jet precession just because of the missing of several 
cycles between adjacent components. However, even under this assumption the acceleration 
of jet precession from the ghost bubbles to the active jets is still significant. Because the 
precession is steady when the source activity changes dramatically, the mechanism for the 
precession is independent of the accretion and thus most likely of the geodetic precession 
or orbital motion of SMBHB. Under the assumption of steady jet precession, we discussed 
the two possible mechanisms. Our discussions suggest that if the precession is due to the 
binary orbital motion, the SMBHB should have a too small mass ratio (g ^ 2 x 10^^) 
and the acceleration of jet precession from the ghost bubbles through outer lobes to the 
active jets cannot be explained reasonably. Our results show that the steady precession 
from the ancient bubbles to the outer lobes is probably due to the geodetic precession 
and the jet precession from the outer lobes to active jets may be due to the precession 
of a warped circumbinary light standard thin disk. In this scenario, SMBHB formed in a 
major merger with mass ratio q = 0.21021^^^^ and the binary has a separation 
a ~ 9.1 X lO^rGCtZi^^ fo^^ S]^2^ — 0.29 pc. The predicated acceleration of jet precession 
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is about a few to ten percent and may have not yet been observed because of the low 
observational accuracy. 

Like what our theoretical investigations show that there is no a priori requirement for a 
steady jet precession, we discussed the implications that if a continuous rapid acceleration of 
precession from the ancient bubbles to the active jets has indeed been observed. Our results 
show that in this case the mechanism for jet precession is the orbital motion and that the 
rapid acceleration of jet precession is due to the rapid evolution of SMBHB because of galactic 
dynamic friction. The calculations give a galactic dynamic friction evolution timescale Ta ~ 
(6 — 9) X 10'' yr and an averaged dynamic friction velocity ^ ^ —1.54 x 10®cm/s. The 
SMBHB forms in the major galaxy merger with an averaged black hole mass ratio q ~ 0.76 
and has a separation a fa 0.8 — 1.46 Kpc. 

As our conclusions, we discussed the scenarios for jet precession in AGNs and calculated 
the time derivatives of the precession timescale. Our calculations show that jet precession 
is accelerated in SMBHB models but nearly steady in the Bardeen-Peterson effect scenario. 
We analytically computed the predicated acceleration of jet precession in the evolution of 
a SMBHB from unbound to gravitational-wave-dominated stages and showed that the time 
variation is significant and can be detected easily. One can estimate the evolution timescale 
and mass ratio of SMBHB, and the parameters of accretion disk in AGNs by measuring 
the central black hole mass, the accretion rate, the jet precession timescale, and its time 
derivative. If we have observations of jet precession acceleration of a sample of radio sources, 
we can test the hierarchical galaxy formation model and the galactic dynamics. We can also 
estimate the fraction of SMBHB that can get coalesced quickly and give rise to gravitational 
wave radiation bursts. 

We are grateful to D.N.G. Lin, J. Magorrian, J.F. Lu, and X.-B. Wu for helpful discus- 
sions and comments. Many thanks are due to the referee for constructive comments, which 
made us to improve the presentation of the paper significantly. This work is supported by 
the National Natural Science Foundation of China (No. 10573001). 
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